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Abstract. Recent experimental data support the presence of quark coalescence 
in heavy ion collisions at RHIC energies. Hadronization of quark matter and 
hadron formation in heavy ion collisions can be described by the coalescence 
process, and measured data are reproduced successfully. On the other hand, the 
theoretical coalescence calculations are based on a non-rclativistic description. 
Here we investigate the robustness of the coalescence description, using different 
wave-function overlap during hadron formation. 

Quark coalescence has been proposed many years ago to describe quark matter 
hadronization [TJ [H [3] • The basic idea was to consider massive constituent quarks 
in the deconfmed phase, which are ready to hadronize through "coalescence", which 
is a clustering process driven by an attractive force between the properly coloured 
quark degrees of freedom. The presence of these massive excitations in the quark- 
matter phase is supported by the analysis of lattice QCD data and the recognition of 
massive quasi-particles in the strongly interacting deconhned phase close to the quark- 
hadron phase transition 4J. The attractive force generated by the gluons (which are 
considered in this way) is modelled by non-relativistic colour potential between the 
quarks and antiquarks. Mesons are produced by quark-antiquark coalescence. Baryons 
are produced in two steps: at first 3 diquarks appear through quark-quark coalescence, 
which is followed by a diquark-quark coalescence into a colourless baryon. 

Although particle yields, ratios and spectra have been reproduced successfully in 
ALCOR [3[6] and MICOR 7, 8 coalescence models, but thermal models were similarly 
successful in the low-p^ region and more widely used because of their simplicity. 
At RHIC energy intense data collection has been performed in the intermediate-pT 
region (3 < px < 8 GeV/c) and the measured anomalous proton/pion ratio could 
have been explained by quark coalescence and recombination models [TOl [TT] . This 
fact increased the interest in this microscopical hadronization mechanism and more 
applications appeared. The recognition of valence quark number scaling in the data on 
asymmetric flow (1)2) supported very strongly the quark matter formation and quark 
coalescence at RHIC and SPS energies [T2] , 

The success of the quark coalescence model raised an interesting question, namely 
why this model is working successfully, when it is based on quantum mechanics and 
non-relativistic quark-quark interaction picture. The main reason is the following: 
although the quarks are moving with a relativistic velocity out of the reaction volume, 
but they can participate in the coalescence process only if their relative velocity is 
small. In this case quantum mechanics becomes valid in the comoving system and 
wave functions start to play a significant role in the description of hadron formation. 
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In this paper we investigate the sensitivity of the hadron yields on the use of the 
different quark and hadron wave functions. We will consider plain wave and gaussian 
wave functions for quarks and antiquarks, where the gaussian choice indicates the 
presence of a finite volume deconfined quark matter. For mesons we use the same two 
choices and extend the list with the hydrogen-like wave function, which is connected to 
an analogy between electromagnetic and strong forces with proper coupling constant. 
For diquarks we use the same three choices. However, the diquark and the plasma 
state may interact and the diquark wave function can be modified, flipping between 
the plain wave, gaussian and hydrogen like wave functions. Thus baryons are produced 
similarly to mesons, but in two steps, as we will discuss it in details. 

At first we determine the quark coalescence cross section in quantum mechanics. 
We use the quantum mechanical pick-up reaction [15] : qi + Q\ q2 — ► h + Q', where 
quark qi picks up quark qi from the plasma Q, producing prehadron h and plasma 
Q' . The quantum mechanical amplitude of the pick up reaction is given: 

9gh = V g ~^ h ' Q ' f d 3 f!d 3 f 2 ^*{2 u S 2 )V{Si ~ x 2 )Mxi)Mx2) , (1) 



2tt 

where (f>i(xi) is the wave function of the qi quark, and ^{x\,X2) belongs to the 
prehadron with mass Mh- Mh.Q> is the reduced mass of h and Q'. Since Mh <C Mqi 
therefore Mh,Q' — Mh- The standard two-body coordinates can be introduced as 
relative distance vector (r), relative momentum vector (k), space and momentum 
vector of the center of mass (X, P). Since the outgoing prehadron has the momentum 
P = pi + P2, thus the wave function "J is simplified to ^(r) • e tPX . 

In eq.([T]) V{x\ — x 2 ) denotes the quark-quark interaction. Here we introduce the 
Yukawa-potential, which depends on the relative distance, r: 

F(r-) = -a(A i A j >^ (2) 

The screening mass has a relatively large value at T as T c : m g = 0.8 GeV [I]. 
The colour factor (AjAj) is determined by the colour combination of the interacting 
particles. In the limit of m g = Coulomb potential is restored, which has been 
used in the ALCOR model 1 . In the MICOR model a Yukawa potential has been 
considered [7]. 

In 2 — > 2 reactions (e.g. a + b — > c + d) the cross section can be determined from 
the amplitude as 

a(k) = ^ \g gh (k)\ 2 , (3) 

Vab 

where v a b and v cc i are the relative velocities. In our case the factor v cc i could be large, 
because it is the relative velocity of the outgoing (pre)hadron and the plasma. Thus 
we have two choices during evaluation of eq. ([3]): 

P 2 

O"0Ki 92 1 classical = \9gh{k)\ (4) 

Mh 

P 2 

o-{k)v qiq2 \ re i at i msU c = , \9gh{k)\ (5) 

V p + M l 

Now we are ready to introduce quark and prehadron wave functions to calculate 
numerically applicable coalescence cross sections from eqs. (HJ-©. In the lack of 
precise quark wave functions, we will use simple functions, which will be displayed 
and discussed after introducing the calculation of hadron production rates and yields. 
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Previously, the coalescence cross sections were obtained from quantum mechanics. 
The wanted hadron production rates and yields can be determined by a statistical 
method based on rate equations using the former cross sections. 

Prehadron h is composed from quarks q\ and q2 via coalescence, and its 
production is proportional to the densities of the constituents, n\ and ^2 [U EI : 

d^{n h u^) = (CT12U12) n 1 n 2 ■ (6) 

The rate, (o^^), ^ s calculated as a phase space average: 

, h y x _ I d 3 Pld 3 p2d 3 X 1 d 3 X2Pl2(x 1 ,X2)fl(xi,pi)f 2 {x2,P2)o-Vl2 ,^ 

12 12 J d 3 p 1 d 3 p 2 d 3 x 1 d 3 X2Pi2(xi,X2)fi(x 1 ,p 1 )f2(x2,P2) 

where fi(x*i,Pi) are the quark momentum distributions and p describes the locality 
of the quark coalescence. Requiring that quarks with the same location are able to 
coalesce, p becomes a simple Dirac delta. Assuming isotrop plasma state, eq.([7]) is 
simplified into the following expression: 

, h v Jd 3 p!d 3 p2 f q (m 1 ,p 1 )f q (m 2 ,P2)avi2 , - 

{ai2Vl2> fd 3 p!d 3 p2 /g(m 1 ,pi)/q(m 2! p 2 ) { ) 

During the evaluation of the rate one can use any proper distribution function. 
Because of the massive quarks we can use non-relativistic Boltzmann distribution: 

ff oltzmann (m,p) = e-^r (9) 

In parallel, relativistic quark distribution can be used also, namely Jiittner 
distribution, which is simplified in the local rest frame of the expanding fireball: 

f^ uttner (m,p) = = P - V " 2 r m2 (10) 

We will use both distributions and investigate the sensitivity of the rate on this choice 
and the temperature T. Furthermore, in the above quark momentum distribution 
functions we will use quark masses m q — 0.3 GeV and m s — 0.5 GeV, which values are 
verified in the analysis of the lattice data close to the quark-hadron phase transition [J] . 

Now we have all expression to investigate prehadron production from quark 
matter. If we define the necessary wave functions for quarks and prehadrons, then we 
can perform numerical calculations and determine the wanted particle productions. 
As we mentioned earlier, we introduce two types of wave functions for quarks: 

<j>i = — ^— e ipiSi PW : plain wave (11) 



1 

I^Tgip.z, G : gaussian (12) 



ri (27rp 2 ) 3 / 4 
For prehadrons we introduce the following wave functions: 

* = — L^Le 1 ^ PW: plain wave (13) 

* = — )=— = e~« e" 3 * H : hydrogen like (14) 

y/Vf V ira 3 

* = G: gaussian (15) 



These expressions contain the volume normalization factors related to the 
characteristic volume of quarks (V q ), prehadrons (Vh), and the fireball volume (V/). 
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Figure 1. The dependence of the coalescence rate ((crv)) on gluon mass (m g ) at 
fix temperature (T = 180 MeV) (left) and on temperature (T) at fix gluon mass 
(m g = 800 MeV) (right) for the case of quark plain waves coalesce into mesonic 
plain wave. 



The wave functions in eqs. (fTTjl - (fT2|) and (fT3 |) -(fl~5 f offer us 9 combinations, which 
can be doubled by the application of Boltzmann and Jiittner momentum distributions. 
For consistency we will calculate meson and baryon productions in the same (fixed) 
wave function combination. We note that ALCOR model pQ used the PW + PW — > H 
combination, M1COR model [7 applied the PW + G — > H choice, including a gaussian 
localization for the picked-up quark. 

In the case of process PW + PW — ► PW averaging with the Boltzmann 
distribution (named as case 'Bl'), one can obtain a compact expression for the rate: 



h \ VgVt M h (m 1 +m 2 ) 2 a 2 s 1 

V'Vh {m 1 m 2 ) a '- i 7T T J (k 1 + m g z y 



This expression shows that the volume terms can be collected into a prefactor together 
with the coupling constant a s , which is valid for every combinations. This prefactor 
should be htted from one data point, and all other particle yields become calculable. 
On the other hand, this factor will drop out from particle ratios (calculated with the 
same wave function combination). In the remaining 17 combinations much longer 
expressions appear for the rates, but all of them can be calculated numerically. 

In our investigation we recognized that combinations containing at least one quark 
plain wave resulted numerically different values, but they differed in a constant factor, 
only. This is the reason why ALCOR's and MICOR's results on particle yields are 
generally the same, although different wave function sets have been used. 

Now we investigate the sensitivity of the coalescence rate on the mass parameter 
of the Yukawa potential (m g ) and the temperature (T). In Figure 1 we display the 
calculated rates (with an arbitrary scale) in the case of PW + PW — > PW, using 
the expression of eq. p^|) . One can see (left) that the meson production rates 
strongly depend on the gluon mass at fix temperature (here we choose T = 180 
MeV). Considering the temperature dependence (right), the rates are very much 
insensitive on this parameter at fix gluon mass (m g = 800 MeV) . This fact verifies the 
applicability of coalescence models for cross-over phase transitions, where no sharp 
transition temperature exists. In Figure 1 one can recognize the presence of a flavour 
dependence, which is connected to the heavier mass of the strangeness. 
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Figure 2. The gluon mass (m g ) dependence of the production ratios of mesons 
with different flavour combinations, relative to the light meson at fix temperature 
(T = 180 MeV). We display the numerical results for the PW + PW -> H (left) 
and G + G -> G (right). 



This flavour dependence is investigated further through the particle ratios, in 
which case uncertainties connected to unknown volume factors disappear. The 
temperature dependence of the ratios is negligible, as one could expect from Figure 1. 
We have found that the strong gluon mass dependence drops out for ratios in all 18 
wave function combinations. In Figure 2 we show our numerical results for the case 
of PW + PW — > H (left) and the case of G + G -> G (right), which results illustrate 
the weak dependence of the ratios on the gluon mass in very different cases. 

These results indicate and prove the robustness of the coalescence model within 
a fixed wave function combination, since the theoretical results depend very weakly 
on temperature and/or gluon mass. This feature nicely support the applicability of 
coalescence models to describe cross-over phase transition, which is expected to happen 
in a wide temperature region. 

Now we would like to investigate the sensitivity of the coalescence model on 
different wave function setups. For this task we want to use real data measured at 
RHIC energy. Our method is the following: we use one measured meson ratio to fix 
the open parameters, recalculate other measurable hadron ratios, and investigate the 
difference between the obtained results in difference wave function setups. One of our 
candidate for the starting point is the ratio §/K* — 0.60 ± 0.15 measured at RHIC 
in central Au+Au collisions at yfs = 200 AGeV [13]. (We use the middle value and 
neglect the influence of the error.) In the coalescence model this ratio appears in the 
following way: 

A = (a h v)ss-N 8 N a 

Nk* (o- h v) sq ■ N s N q [ ' 

The rates in eq. (|17p can be calculated in a fixed wave function setup. Thus we can 
determine the missing N s /N q factor from the measured N^/Nk* value. Then different 
strange-non-strange meson and baryon ratios can be calculated from the model. Here 
we will calculate the ratios K* /p°, E*/A, S /£ and f2/H . These resonances have been 
measured [15] and the above ratios could be determined from existing experimental 
data. 
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Table 1. Hadron ratios in different wave function setups with Boltzmann (upper 
part) and Jiittner distributions (lower part). The ratio N s /N q is determined from 
the fixed ratio Q/K* = 0.6 measured at RHIC | 14| . Other particle ratios are 
calculated from the corresponding wave function setup, using the former strange 
to light quark ratio. 



Model 




iV 3 /iV 8 


K*/p° 


S*/A 




n/s* 


PW+PW -» 


PW 


0.448 


0.671 


0.790 


0.700 


0.644 


PW+PW -» 


H 


0.235 


0.598 


0.848 


0.785 


0.746 


PW+PW -> 


G 


0.474 


0.686 


0.791 


0.690 


0.628 


G+G -> PW 




0.441 


0.668 


0.790 


0.702 


0.648 


G+G -> H 




0.227 


0.593 


0.849 


0.791 


0.753 


G+G -> G 




0.503 


0.700 


0.791 


0.681 


0.614 


PW+PW 


PW 


0.401 


0.701 


0.857 


0.721 


0.644 


PW+PW -♦ 


H 


0.206 


0.624 


0.928 


0.815 


0.752 


PW+PW -> 


G 


0.424 


0.706 


0.843 


0.708 


0.631 


G+G -> PW 




0.396 


0.701 


0.862 


0.724 


0.647 


G+G -> H 




0.201 


0.620 


0.931 


0.819 


0.757 


G+G -> G 




0.435 


0.712 


0.845 


0.706 


0.627 



Table 1. summarizes our numerical results. We can see that different wave 
function setups result very different N s /N q values, namely 100-150 % difference can 
be seen between the smallest and the largest values of quark ratios. On the other 
hand, this uncertainty drops to a 10-15 % difference, both for mesonic and baryonic 
ratios. These results prove most strikingly why the coalescence models yield very good 
agreement during data reconstruction, if we start from one measured values. Further 
analyses are in progress to reveal a deeper connection between the manifestation of 
conservation laws and the structure of the quantummechanical description of quark 
coalescence. 
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